Exact solution of Z 2 Chern-Simons model on a triangular lattice 
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We construct the Hamiltonian description of the Chern-Simons theory with Z n gauge group on 
a triangular lattice. We show that the Z2 model can be mapped onto free Majorana fermions and 
compute the excitation spectrum. In the bulk the spectrum turns out to be gapless but acquires 
a gap if a magnetic term is added to the Hamiltonian. On a lattice edge one gets additional non- 
gauge invariant (matter) gapless degrees of freedom whose number grows linearly with the edge 
length. Therefore, a small hole in the lattice plays the role of a charged particle characterized by a 
non-trivial projective representation of the gauge group, while a long edge provides a decoherence 
mechanism for the fluxes. We discuss briefly the implications for the implementations of protected 
qubits. 



I. INTRODUCTION. 



The challenge of the error free quantum computa- 
tio n 1 ! 2 ! 3 resulted in a surge of interest to many physical 
systems and mathematical models that were considered 
very exotic before. While it is clearly very difficult (if not 
impossible) to satisfy the conditions of long decoherence 
rate and scalability in simple physical systems^, both 
can be in principle satisfied if elementary bits are rep- 
resented by anyons, the particles that indergo non-trivial 
transformations when moved adiabatically around each 
other (braided)£*§*£. One of the most famous examples of 
such excitations is provided by the fractional Quantum 
Hall Effects^. The difficult part is, of course, to iden- 
tify a realistic physical system that has such excitations 
and allows their manipulations. This problem should be 
separated into different layers. The bottom layer is the 
physical system itself, the second is the theoretical model 
that identifies the low energy processes, the third is the 
mathematical model that starts with the most relevant 
low energy degrees of freedom and gives the properties of 
anyons while the fourth deals with construction of the set 
of rules on how to move the anyons in order to achieve a 
set of universal quantum gates (further lies the layer of 
algorithms and so on). 

One of the most interesting set of problems of the third 
layer is provided by the Chern-Simons theories with dis- 
crete groups on the lattice: on one hand in these theories 
one expects to have a non-local interaction between fluxes 
that gives them the anyonic properties, on the other hand 
they might describe physics of some solid state arrays. 
In particular, we have shown very recentljiiS that TL-i 
Chern-Simons theory can be realistically implemented in 
a Josephson junction array on a square lattice. Further- 
more, the ground state of these physical systems is dou- 
bly degenerate but locally instinguishable so that even a 
small sized lattice provides a very good protection of this 
degeneracy against the external noise. Such protection is 
expected to become ideal for the larger lattice sizes if the 
gap to excitations is finite. The lack of the exact solu- 



tion did not allow us to make definite conclusions on the 
properties of this model for larger lattice sizes. However, 
the analytical solution in limiting cases combined with 
the extensive numerics^ indicates that the gap for the 
low energy excitations closes in the thermodynamic limit. 
Practically, this would imply that protection that can be 
achieved in such arrays does not grow with the array size 
beyond a certain limit. In order to get a better under- 
standing of the properties of such models we consider 
here a similar (2 n Chern-Simons) model on a triangular 
lattice. It turns out that for n — 2 the model can be 
solved exactly by a mapping to Majorana fermions; we 
find the gapless spectrum of the photons which acquires 
a gap if additional ('magnetic') terms are included in the 
Hamiltonian. We hope that long wave properties of this 
model are analogous to the properties of the square lat- 
tice model and thus the problem of the low energy exci- 
tations of the latter can be remedied by the 'magnetic' 
field term in the Hamiltonian that plays the role of a tun- 
able chemical potential for fluxes. Because it is difficult 
to avoid boundaries in a physical implementations, it is 
important to understand what is the effect of edges on 
the Chern-Simons theory. In particular, from the view 
point of the topological protection, it is important to un- 
derstand whether they lead to gapless excitations such as 
edge states in Quantum Hall Effect. 

Generally, while the lattice gauge theories without 
Chern-Simons term are well studied and understood, 
much less is known about their Chern-Simons counter- 
parts. The reason for this is that Chern-Simons term 
implies coupling of the charge and flux; on a lattice the 
charge resides on the sites while the flux is associated 
with the plaquettes. Furthermore, one usually wants to 
preserve the symmetry of the lattice. To satisfy these 
criteria one has to couple charge with the total flux of 
a few adjacent plaquettes which leads to novel features 
absent in the continuous theories. We discuss the de- 
tails of this construction in Section II. In particular, for a 
continuous Abelian gauge group the Chern-Simons term 
remains quadratic but becomes non-local in space, i.e. 
its Fourier transform acquires momentum dependence; 
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moreover it is zero for some values of the momenta in 
the Brillouin zone. This leads to the appearance of the 
gapless modes in these theories in the absence of magnetic 
energy. Our analysis presented in Section III of the ex- 
actly solvable Chcrn-Simons theory with TL2 group shows 
that this qualitative picture holds for the discrete group 
as well. Finally, in Section IV we study the boundary 
effects and find that they are also similar for continuous 
and discrete groups. Namely, for continuous groups, the 
general arguments show^ the appearance of a non-gauge 
invariant matter field on the boundary in Chern-Simons 
theories, similarly the exact solution of 7Li model shows 
that the same phenomena happens in discrete models as 
well. 



II. MODEL WITH Z n SYMMETRY 

Let us first begin to discuss the construction of the 
Chcrn-Simons model with Z„ symmetry on the triangu- 
lar lattice. The main element of this construction are 
the expressions for the electric field operators, and 
the local field translation operators . 10 Both should 
preserve the symmetry of the lattice; furthermore, the 
expression for the electric field operators, £^ should be 
gauge invariant while those for should preserve the 
electric fields operators. By analogy with our previous 
discussion for the square lattice^, we write: 
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Here Am and II mn are (discrete) gauge potential and 
canonical conjugate momentum on lattice link that sat- 
isfy usual canonical commutation relations: 
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where (kl) and (mn) denote oriented links on the trian- 
gular lattice. Reversing the orientation of a link changes 
the sign of the corresponding variables A and II. In the 
case of a 2„ symmetry, the local vector potentials A k i 
are constrained to be integer multiples of 2-njn. The 
Chcrn-Simons coefficient v(kl] mn) is zero unless (mn) 
is one of the four neighbor links of (kl) as illustrated on 
Fig-H In this case, v(kl\ mn) = if the link (mn) is ori- 
ented from right to left for an observer standing on link 
(kl) and looking towards site I. For the converse relative 
orientation, v(kl\ mn) = —v 

The main consequence of these definitions is that 
the commutation relations of electric field operators on 
nearby links are modified by phase factors. For (kl) and 
(mn) oriented as on Fig. ^ we have: 
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FIG. 1: (Color online) The orientation convention used in the 
construction of a Chern-Simons model on a triangular lattice. 
The thick (black) line represents a Tiki operator, conjugated 
to the vector potential attached to the link. The arrow is 
required since Tiki = — Iljijfe. The gray (red) lines keep track 
of signs of u(kl; Im) coefficients which appear in the phase- 
factors imposed by the Chern-Simons term in the definition of 
local electric operators £^ and field translation operators . 
Gray (red) lines are oriented from site m to site n whenever 
the coefficient v(kl; mn) is positive. 



Similarly, the local field translation operators satisfy: 

,2tt 



T k JT+ n = exp -i(—Y2v T+ n T k 
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We shall focuss here on the case where these phase factors 
are equal to — 1. This requires v to be of the form: 



, = n(-y(m + -) 



(3) 



where m is an integer. 




FIG. 2: (Color online) Construction of the local gauge trans- 
formation generator based at the central site, see Eq. @. 
Thick full links correspond to electric field operators E^, and 
grey (red) oriented links to phase factors, exp(i^^-A mn ), 
whose total contribution is equal to exp(— i^^<3>fc), where 
<E>fc is the total flux (counted counterclockwise) through the 
hexagon. 

To construct the Z n model, we start from the sys- 
tem where local vector potentials can be arbitrary (and 
in particular unbounded) integer multiples of 27r/n, and 
keep then only the subspace of periodic states which 
satisfy: 
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In general, this condition is not compatible with the pos- 
sibility to perform arbitrary local gauge transformations, 
because the T kl no longer commute with generators of 
local gauge transformations Uk defined as: 



(fc) 



(fc) 



(5) 



where §k = Tlhex Amn is the total flux through the loop 
composed by the six first neighbors of site fc and ori- 
ented counterclockwise, as illustrated on Fig. This is 
an integer multiple of 2n/n. Note that some care is re- 
quired in chosing the order of the six operators involved 
in the above products. A convenient choice is to lump 
together each link (kl) with its opposite (hi 1 ). This con- 
vention is sufficient to specify the total phase factor be- 
cause 71j operators on all but adjacent links commute (so 
that [T$, T£,\ = 0) while different pairs Thi^kl' commute 
because of the cancellation of the phase factors for any 
value of v. So, this choice yields three pairwise products 
which commute among themselves. With the notations 
of Fig. |2 we have: 



Uk%n 



exp 



) T mn U k ((») 
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Therefore, the 2„ periodicity conditions are compati- 
ble with local gauge invariance provided the quantity 
(27r/n) 2 4^7i is an integer multiple of 2ir, which is clearly 
the case for v chosen according to Eq. ©. Another im- 
portant consequence of this choice is the very simple con- 
nection between generators of gauge transformations and 
local electric fields: 



(fc) 



U k =l[£ 



A last consequence is that for any even value of n: 
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Once this set of basic operators has been defined, the 
natural gauge-invariant Hamiltonian is constructed as a 
sum of two terms: the first one involves local electric field 
operators on individual links and the second one the local 
magnetic fluxes around plaquettes. Specifically: 



(9) 
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The first sum is over links (k,l) and the second one over 
elementary triangular plaquettes (j, fc, I) which are ori- 
ented counterclockwise. At this stage, / may be any 
function and the flux 4>(j, fc, I) = Aij + Ajk + Am- Note 
that we are not dealing here with the pure Chern-Simons 
theory, but rather with a discrete analogue of a Maxwell- 
Chern Simons theory. Indeed, the former has only a 
small Hilbert space which dimension is independent of 



the system size, and a vanishing Hamiltonian. In the pure 
Chern-Simons theory, only fluxless configurations are al- 
lowed, unless the ambiant space has a non-trivial topol- 
ogy. The Hilbert space of the Maxwell Chern-Simons 
theory is much larger, and is much more likely to corre- 
spond to the low energy sector of a real physical system. 
The ground-state sector of this model is then expected 
to be described by a pure Chern-Simons model. 



III. THE n = 2 CASE 

From now on, we assume n — 2. We shall map the 
subspace of gauge invariant and 2tt periodic states (the 
"physical subspace") on a system of Majorana fermions 
attached to the plaquettes of the original triangular lat- 
tice. The periodicity condition implies, thanks to eq. @ 
that physical states \^f>) satisfy: 



(^) 2 |*> = l*> 



We may then set: 
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S kl = £ k i 



(10) 
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on this physical subspace. As for the square lattice, it is 




FIG. 3: (Color online) Strings of electric field operators (thick 
lines), S r and S T > that create or remove IL2 fluxes in the pla- 
quettes indicated by the shaded triangles. For two arbitrary 
fluxes the corresponding strings have a common part (full 
lines) and individual parts (gray lines) . Two individual parts 
form a contour that connects two fluxes. Each of these three 
substrings anticommutes with all others that results in the 
anticommutation of the flux creation operators (see text) 

fruitful to introduce string-like operators S r which create 
or remove a W.2 flux on a plaquette labelled by a site r 
on the hexagonal dual lattice. This operator involves the 
product of Ski on a contour that starts from a specific 
bond of the lattice, as shown on Fig. |SJ For the gauge 
invariant states all such contours are equivalent, so this 
contour can be chosen arbitrarily. For instance, one can 
choose the contours that first go up and then left toward 
the plaquette r. Using the fact that electric operators on 
nearby links anti-commute, we have: 



5 r 2 = (-1) 



l(r)-l 



(12) 
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where l(r) is the length of the contour r, namely it is 
the total number of electric operators involved in the 
construction of S r . Because a gauge transformation 
changes the length of the contour by an even number, 
the right-hand side of Eq. I|12l) depends only on the par- 
ity [r] = (l(r) — l)mod2, of the plaquette r and not on 
the contour leading to it. Note that S r is not always a 
hermitian operator, since 



Si = {-l) [r] S r 



(13) 



The commutation rules obeyed by these string operators 
are the following: 



{S r ,S r ,} = 26(r,r')(-l)^ 



(14) 



This crucial property which exhibits a clear fermionic be- 
havior can be proved in two ways. First, one can deform 
the contours so that they overlap between point and r. 
Then S T ' — S r ^ r S r where S r ^ r is the contour beginning 
at r and ending at r'. This contour S r i_ r contains exactly 
one electric field operator namely its first one, which anti- 
commutes with the last electric field operator entering in 
S r , so operators S r ^ r and S r anticommute which proves 
(|14H . Alternatively, one can write the original S opera- 
tors as products S r — S r So and S r > = S r >So where Sq 
is the product of electric field operators on the common 
part of the contour leading to plaqucttcs r and r', see 
Fig- El The operators S r and So have exactly one pair of 
nearest neighboring electric fields, so they anticommute: 
{S r ,,S } = {S r ,S } = 0. Thus, S r S T/ = -S$S r S r , and 
S r >S T = —SySr/Sr, further, noticing that the operators 
S r S r > contain exactly one pair of anticommuting electric 
fields, we get l|T4"|). 

A very important property of this model restricted to 
its physical subspace defined above is that local electric 
operators on link (kl) are simply related to bilinear ex- 
pressions of the form S T S T ' where the plaquettes r and 
r' are located on both sides of the link (kl). If r and r' 
are nearest neighbor plaquettes, we have: 



(15) 



which is a direct consequence of the definition of string 
operators. Thus, 



Ski = (-i) [r] s r s r 



(16) 



The last stage is to introduce Majorana fermions Xr 
related to string operators by: 

Sr=l [r] Xr (17) 

With this definition, it is easy to check that: 



Eq. becomes: 



Xr = Xr 
{Xr,Xr>} = 25 r 



£kl = IXrXv 



(18) 
(19) 



(20) 



if site r' is even ([r'] = 0) and site r is therefore odd. 

The electrical part of the Hamiltonian thus maps into 
the hopping of Majorana fermions. The full Hamilto- 
nian includes also the magnetic field part. In case of Z2 
model it is given by the operator that takes two values 
depending on the flux in a given plaquette, so a general 
function /($) from @ is reduced to a linear function 
/($) = Writing the Majorana fermion as a sum of 
two usual fermion operators x = c + c ) we see that the 
total fermion number e'e changes by ±1 by the flux cre- 
ation operator, so in the fermion language the full Hamil- 
tonian becomes: 



H 



1 Yl ( C r + 4)( c s + 4) + V>y] 4°r 



(r,s 



where the indices r, s run over the sites of the dual 
(hexagon) lattice and the first sum goes over all nearest 
neighbours on this lattice. The spectrum of the excita- 
tions 



E(k) = Vv 2 + \t(k)\*±\t(k)\ 



(21) 



where t(k) is spectrum of fermions on the honeycomb 
lattice with a purely nearest neighbour hopping Hamil- 
tonian H t — J2( r s ) c it c s- To compute t(k) we choose ele- 
mentary cell consisting of two sites, for instance the ones 
that belong to the same vertical bond. In momentum 
space the Hamiltonian becomes a matrix H = c k K a bCb } k 
with 



K 







1 + e ik « + e ik?? 



where £ = (-^, §), 77 = (— |) are unit vectors of the 
lattice of vertical bonds. We get 



((*) 



\ 



1 + 4 cos 





In the absence of magnetic term the spectrum has 
dispersionless mode and a dispersive one which has 
zero only at the isolated Fermi points (±4tt/(3-\/3), 0), 
(±27r/(3v3), 27r/3). Near the Fermi point the spectrum 
is linear E(k) = 3\k\. In the presence of magnetic term 
two things happen: the dispersive mode acquires a gap 
2/ii, near the Fermi point the spectrum becomes massive 
E(k) = v /(3/2) 2 |fc| 2 +7t+(3/2)|fc| and the dispersionless 
mode acquires /c-dependence with the gap \/ [i 2 + 9 — 3. 

The presence of a dispersionless band for /1 = is di- 
rectly connected to the presence of a set of local symme- 
tries in the model. Indeed, if we introduce the Majorana 
operators x = i(c— c'), they satisfy: 

(22) 
(23) 
(24) 



Xl = Xr 



{Xr,Xr'} = 2<S r , r 
{Xr,Xr} = 



As a result, Xr commutes with H for any r, therefore 
providing a set of non-commuting local symmetries which 
are destroyed by the presence of a magnetic term ji. 
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IV. EDGE STATES 



An interesting feature of the present model is that it 
exhibits edge states localized around boundaries. This is 
known to be a general property of Chern-Simons theo- 
ries^, but an advantage of a lattice gauge theory with 
a discrete gauge group over a model defined on continu- 
ous space is to yield a finite dimensional Hilbert space. 
Therefore, we do not require any sophisticated analy- 
sis to specify boundary conditions. In the presence of a 
boundary, the definition of basic operators and 
for a link kl along the boundary is the same as before, 
with the exception that the associated phase-factor only 
involves a smaller set of vector potentials A mn attached 
to the neighboring links lying inside the system. It is 
then easy to check that this does not modify the basic 
commutation relations between these operators. In par- 
ticular, we shall still concentrate here on the n = 2 case 
where nearby local field translation operators anticom- 
mute. The gauge generators Uk located at site k is still 
the product of all Tir operators such that kl lies inside 
the system. As illustrated on Fig. 0] this implies that 
Uk and Ui still commute when at least site k or I is not 
on the boundary, but they now anticommute if k and I 
are nearest neighbor sites both located on the boundary. 
This holds for most possible shapes of the boundary, with 
few exceptions which are depicted on Fig-EI The striking 
consequence of this is that it is no longer possible to di- 
agonalize simultanenously all local gauge generators Uk 
belonging to the same boundary. In particular, gauge sin- 
glets no longer exist since they are replaced by degenerate 
multiplets corresponding to projective representations of 
the gauge group associated to boundary sites. 




180 60 



FIG. 4: (Color online) Commutation relations between gauge 
generators attached to boundary sites in the Chern- 
Simons model. We show various possible edge shapes and in- 
dicate the corresponding opening angle expressed in degrees. 
For each situation, we display the gauge generators for a pair 
of nearest-neighbor sites. These generators involve the elec- 
tric operators originating from a given site, and are depicted 
either as full (red) thick lines or dashed (blue) thick lines. 
For all situations appearing on the outer boundary of the lat- 
tice, these gauge generators are found to anticommute, which 
induces degenerate multiplets for the gauge symmetry. An 
exception to this rule is the —60 cusp shown on an inner hole 
for which the two nearby generators commute. 



numbers corresponding to: 

U 2 n\V) =r 2n |tf) (25) 

L 

I] =r|tf) (26) 

71=1 



Let us now describe these multiplets in some detail. 
For this, we consider a finite triangular lattice with an 
outer boundary, and possibly with some inner holes, each 
of them bringing its own boundary. We assume that 
these holes are not too close from each other, so that lo- 
cal gauge generators associated to sites belonging to two 
different boundaries always commute. This assumption 
allows us to treat each boundary separately from the oth- 
ers. To simplify the discussion, we shall consider first the 
edge without —60 degrees turns that contains an even 
number L of sites, which will be labelled by indices n 
running from 1 to I. We suppose that fluxes in each 
plaquette and through each inner hole have been fixed, 
thereby concentrating on the degrees of freedom attached 
to gauge transformations. We may then associate to 
our boundary a Hilbert space containing 2 L independent 
states. In order to construct irreducible representations 
of the boundary gauge-group, we have to find a maxi- 
mal subset of mutually commuting generators, which are 
then diagonalized simultaneously. Such set may be cho- 
sen as follows: take first generators U 2n on even sites, 
and add then a global gauge generator I7n=i ^™ f° r ^ ne 
boundary. This yields 2 1+L / 2 possible sets of quantum 



where eigenvalues T2 n and r can be ±1. We may still 
interpret r as the total 7Li charge of matter induced 
on the boundary by the Chern-Simons term. Starting 
from a common eigenstate |\&) as above, applying L^n+i 
produces a new eigenstate in which Tin and T 2n +2 have 
changed sign simultaneously. Note that we have: 

2 2 
IJ C^-il*) =rj]r 2 „|*) (27) 

n—l n—1 

so only (L/2) — 1 generators located on odd sites act inde- 
pendently. We therefore generate a 2^ i / 2 ^ 1 dimensional 
irreducible multiplet, so the Hilbert space attached to the 
boundary splits into 2^ L / 2 ^ +1 such degenerate subspaces. 

The above construction provides a basis in each mul- 
tiplet for which U2 n is represented by a diagonal matrix, 
whereas l/an+i's play the role of raising or lowering oper- 
ators. A convenient way to visualize these multiplets is 
to describe them in terms of an effective spin 1/2 model, 
attached to even boundary sites. This correspondence is 
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given by: 

U 2 n = ^ n , (l<n<§) (28) 

Ui = tIt* ' (29) 

tfan-i = 7?„_ a 7? n , (2<n<~-l) (30) 

^-i = ^fll^„jr£_ 2 r£ (31) 

A given multiplet corresponds to a fixed eigenvalue for 

Iln=i r 2ni so that the number of independent Ising spins 
is only (L/2) — 1 as discussed above. Note that these 
conclusions apply to sufficiently large holes in the lat- 
tice. For very small holes consisting of two triangles 
glued together to form a rhombus (and of course for a 
single triangle itself) all boundary operators commute. 
The first pair of anticommuting operators appears in a 
trapezoidal hole shown in Fig. This hole carries an 
effective spin 1/2 degree of freedom because it is possi- 
ble to diagonalize all gauge generators except one. Note 
that in this special case the product Jln=i U n anticom- 
mutes with the two local generators on upper sites, so the 
corresponding charge is no longer a conserved quantum 
number. In larger holes, such as hexagon, all nearest 
neighbor gauge generators anticommute. For instance, 
the boundary matter of elementary hexagon can be de- 
scribed as an effective system composed of two spins 1/2. 

What is the interaction between these new matter-like 
degrees of freedom and the fiuxons described in the pre- 
vious section? Suppose first that these fiuxons are not 
allowed to jump accross boundaries. This corresponds 
to a Hamiltonian where the electrical operators attached 
to boundary links are removed. The only interaction 
between the two sets of degrees of freedom occurs via 
the usual Aharonov-Casher phase-factor for fiuxons going 
around inner holes. Each hole generates an effective ~2-2 
orbital magnetic field which produces an adiabatic phase 
4> related to the Z2 charge by e 1 ^ = r. This effect does 
not disturb the internal states of the boundaries which 
remain constants of motion. A stronger interaction oc- 
curs when fiuxons are allowed to cross boundaries (that 
is to disappear from a boundary placquette). A simple 
description of this interaction is possible if a fluxon moves 
around a closed loop which crosses twice a given bound- 
ary. In this process, the final internal state of the bound- 
ary is connected to the initial one by applying the prod- 
uct of all local boundary gauge generators which are also 
enclosed in the closed loop. In the example of a hole with 
the shape of an elementary hexagon, these scattering pro- 
cesses can be simply expressed as products of operators 
taken in the set rf , rf , rf , T4 . In connection with the idea 
of topological quantum computation, we may wonder if 
these well defined operations can be used as a physical 
basis for qubits. The problem with this idea is that lo- 
cal gauge symmetry operators induce transitions between 
the various states of this effective spin system attached 



to the boundary. Since in any real implementation, noise 
acts through basically any local operator, such states can- 
not be protected from environment-induced decoherence. 
If we now consider the effect of either a large boundary, 
or many small inner holes, on the fluxon dynamics, we see 
that these processes tend to entangle fluxon states with 
an effective environment associated to boundary matter, 
resulting in the absence of long-ranged phase coherence 
for the Majorana fermions describing the fiuxons. It is 
also very unlikely that the model preserves its integra- 
biblity, as soon as several fiuxons and boundary matter 
are simultaneously present. 



V. CONCLUSIONS 

In this paper, we have solved exactly a Z2 gauge the- 
ory with a Chern-Simons term on a triangular lattice by 
mapping the gauge-invariant subspace into a free Majo- 
rana fermion system. These fermions keep track of the 
dynamics of local fluxes. In the absence of an energy 
cost to create a Z2 flux, the spectrum is very degener- 
ate: it exhibits a flat band at zero energy and a disper- 
sive one with a vanishing energy at isolated points on 
the Brillouin zone boundary. This picture is very simi- 
lar to the solution found by Kitaev for a quantum spin 
1/2 model on the hexagonal lattice^. In this model, 
he considered nearest-neighbor interactions of the form 
((Xj.7ijj)(<r j-Tiij) where = e x ,e y ,e z depending on the 
unit vector joining sites i and j. These operators have 
the same anticommutation properties as the local fluxon 
moving operators £ki = iXrXr' introduced in the Chern- 
Simons lattice gauge model. However, we have not found 
a way to establish a one to one mapping between these 
two models. We do not see any operator in Kitaev's spin 
model which could be interpreted as a local fluxon num- 
ber, since it would have to anticommute with the three 
components <r x , a v , a z of the local spin at any site. In 
our gauge theory, fiuxons have to be created by non-local 
operators S r associated to contours, so it does not seem 
easy to express them in terms of local spin operators. 

The presence of a gapless spectrum in the absence of a 
magnetic energy term is reminiscent of the gapless phase 
found on a square lattic e) 10 ! 11 . On a square lattice, this 
TLi Chern-Simons model can be mapped into a spin 1/2 
model with anisotropic exchange interactions of the form 
crfcr? or of <r| depending on the orientation of the link 
ij. Adding a magnetic term in the Chern-Simons theory 
is then equivalent to imposing a magnetic field along the 
y direction in the spin model. It would be interesting 
to check that it also induces a spectral gap. Unfortu- 
nately, such a magnetic field term does not commute with 
the non-local symmetry operators associated to rows and 
columns of the square lattice which were responsible for 
the two-fold degeneracy of all energy eigenstates. 

Finally, we have shown explicitely that charged matter 
degrees of freedom are always induced along edges asso- 
ciated to the external boundary or to inner holes. These 
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degrees of freedom appear naturally since the presence 
of an edge forces the representation of the local gauge 
symmetry to become projective, i.e. gauge generators 
attached to nearest neighbor sites both located along an 
edge do not commute. We have shown that these new de- 
grees of freedom provide static Aharonov-Bohm fluxes for 
the orbital motion of fluxons, provided the latter are not 
allowed to cross edges. Transitions within these degen- 
erate matter multiplets are induced by processes where 
a fluxon goes back and forth accross a boundary. Unfor- 
tunately, these multiplets are very sensitive to external 
noise acting through local operators, so it is unlikely they 
may serve the purpose of designing protected qubits. 

It is an immediate extension of this work to add local 
static charges in the system and see how they interact 
with this fluid of Majorana fermions associated to flux- 
ons. Very likely, the interaction will be of the Ahoronv- 
Casher type, which is too weak to induce a confining force 



for a pair of opposite charges. So again, we have a direct 
evidence that a Chern-Simons term distroys the confining 
regime expected in gauge theories with small magnetic 
energy in the absence of a Chern-Simons termM^. 

The most important question connected to protected 
quantum computation is the construction of similar mod- 
els for a finite non-Abelian group££ This is the subject 
of a forthcoming work. 
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